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To clarify the mechanism of recently reported, ambipolar carrier injections into quasi- 
one-dimensional Mott insulators on which field-effect transistors are fabricated, we employ 
the one-dimensional Hubbard model attached to a tight-binding model for source and drain 
electrodes. To take account of the formation of Schottky barriers, we add scalar and vector 
potentials, which satisfy the Poisson equation with boundary values depending on the drain 
voltage, the gate bias, and the work-function difference. The current-voltage characteristics 
are obtained by solving the time-dependent Schrodinger equation in the unrestricted Hartree- 
Fock approximation. Its validity is discussed with the help of the Lanczos method applied 
to small systems. We find generally ambipolar carrier injections in Mott insulators even if 
the work function of the crystal is quite different from that of the electrodes. They result 
from balancing the correlation effect with the barrier effect. For the gate-bias polarity with 
higher Schottky barriers, the correlation effect is weakened accordingly, owing to collective 
transport in the one-dimensional correlated electron systems. 
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1. Introduction 

Molecular materials have been extensively investigated so far. Intriguing properties have 
been observed recently in non-equilibrium conditions. To seek new functions of molecular 
materials that reflect collective characters in non-equilibrium conditions, we here study the 
field-effect properties of one-dimensional Mott insulators. 

Among the electronic devices, the metal-oxide-semiconductor field-effect transistor (MOS- 
FET) is the most important for very-large-scale integrated circuits. Its current-voltage char- 
acteristics are often explained rather simply. One of the conditions that allow the simple 
explanation is the so-called gradual channel approximation, which requires the fact that the 
transverse field (perpendicular to the insulating layer) in the channel is much larger than the 
longitudinal field (parallel to the channel). 1 The total charge per unit area induced in the 
semiconductor is given by the potential difference between the gate electrode and the respec- 
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tive position along the channel, which is multiplied by the gate insulator capacitance per unit 
area. The conductivity is then given by the charge density multiplied by the mobility, which 
is usually regarded as a constant. 

Some field-effect transistors fabricated on organic semiconductors /insulators have very dif- 
ferent situations from the above. For example, carbon nanotube field-effect transistors operate 
as unconventional Schottky barrier transistors, in which transistor action occurs primarily by 
varying the contact resistance rather than the channel conductance. 2-5 Although the gate volt- 
age was initially assumed to modify the nanotube conductance in analogy with the channel of 
an ordinary field-effect transistor (operation involves gate modulation of carrier density in the 
channel), Schottky barriers at the contacts are now known to play a central role (operation 
involves gate modulation of the Schottky barrier height /width) . The Schottky barriers are 
sensitive to the work-function difference between the channel and the source/drain electrodes. 
By matching their work functions (reducing the work function of the electrode by exposing it 
to hydrogen 5 or increasing it by oxygen exposure 2 ), the ambipolar field-effect characteristics 
are achieved. The one-dimensional structure and the coherent transport in carbon nanotubes 
make the contact resistance relatively important, so that the conductivity is no more given 
by the charge density multiplied by the mobility. 

Many molecular materials have strong anisotropy in electronic conduction owing to their 
molecular arrangements in the crystals. In addition, electron correlations also contribute to 
the confinement of the motion of carriers. 6 Recently, Hasegawa et al. reported ambipolar 
field-effect characteristics in metal-insulator-semiconductor field-effect transistor (MISFET) 
device structures based on organic single crystals of the quasi-one-dimensional Mott insulator 
(BEDT-TTF) (F 2 TCNQ) [BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene, F 2 TCNQ=2,5- 
difluorotetracyanoquinodimethane]. 7 Here again, the Schottky barriers play an important role 
owing to the molecular arrangement. In the crystals, the molecular long axes are perpendicular 
to the gate insulator. The direction of the largest intermolecular overlap, which corresponds to 
a transfer integral of about 0.2eV between the neighboring BEDT-TTF molecules, is parallel 
to that of the drain current. 8 The transfer integral between the neighboring BEDT-TTF 
and F 2 TCNQ molecules parallel to the gate insulator is about O.OleV, 8 while the transfer 
integrals along the molecular long axes are much smaller. Thus the motion of carriers is quasi- 
one-dimensional, and especially the motion perpendicular to the gate insulator rarely takes 
place. 

In contrast to the situation involving band insulators, the effect of the interface barrier 
potentials between the Mott insulator and the source/drain electrodes on electron injections 
is found very similar to that on hole injections at the source and drain contacts. 7 Namely, 
the injections of electrons and holes are almost equally possible down to low temperatures. 
The work function of the source/drain electrodes is not needed to match with that of the 
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crystal if it is generally realized. In this paper, we theoretically show that the ambipolar 
field-effect characteristics generally result from combined effects of the electron correlation 
and the Schottky barriers. We calculate the current-voltage characteristics using the one- 
dimensional Hubbard model for Mott insulators (and a tight-binding model with alternating 
transfer integrals, alternating site energies, or staggered magnetic fields for band insulators) 
and adding to it scalar and vector potentials that originate from the long-range Coulomb 
interaction and are modified by the work-function difference, the gate bias, and the drain 
voltage. The time-dependent Schrodinger equation is combined with the Poisson equation 
and numerically solved self-consistently at each site and time within the unrestricted Hartree- 
Fock approximation. 

The ambipolar field-effect characteristics are realized in a wide parameter space spanned 
by the work-function and bandwidth differences between the crystal and the electrodes as well 
as the strength of the long-range Coulomb interaction governing the potential distribution. 
This is shown to be achieved by balancing the correlation effect with the barrier effect. For the 
gate-bias polarity with higher Schottky barriers, the correlation effect is weakened accordingly. 
The charge is transported in a very collective manner in one-dimensional Mott insulators 
in strong-enough electric fields. Such collective properties are reminiscent of the knock-on 
mechanism in ion channels. 9 ' 10 

2. One-Dimensional Model for a Crystal Attached to Electrodes 

Because the carrier motion perpendicular to the gate insulator rarely takes place owing 
to the negligible intermolecular overlaps along the molecular long axes, we assume that the 
effect of the gate electrode appears only in a boundary value of the Poisson equation for the 
source-drain potential. This "boundary" is actually set at the middle of the crystal below. 
Electrons or holes are injected at the source and drain electrodes. Because they move along 
the crystal and these electrodes, we need to include them explicitly in the model. 

We use the one-dimensional Hubbard model for Mott insulators attached to the source 
and drain electrodes represented by a tight-binding model. The total number of electrons is 
the same as the number of sites. In general, the work function of the crystal is different from 
that of the source and drain electrodes. The work-function difference <p is so defined that <fi > 
if the crystal has a higher work function than the electrodes. The model is written as 



where c icr (Q iCr ) creates (annihilates) an electron with spin a at site i, ni a = c icr Ci ta , and 



H 






(1) 
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rii = X^o- n «o-- The site energy ej is set at —<p in the crystal and at in the electrodes. The 
absolute value of the transfer integral | | is set at t c if either i or i + 1 is in the 

crystal and at t e otherwise. The on-site repulsion U{ is set at U in the crystal and at in the 
electrodes. 

Later, we consider band insulators, for comparison, by using a tight-binding model (U=0) 
with alternating transfer integrals [t c above is replaced by i c — (— l) l <5i], alternating site energies 
[—(f) above is replaced by —(f) + (—l)' l 5e\, or staggered magnetic fields [—<f) above is replaced 
by —(f) + (—l) l h for up-spin and —(f) — { — l) l h for down-spin electrons]. 

Hereafter, the drain voltage (relative to the source) is denoted by Vd, and the gate voltage 
(relative to the source) by Vg- Following ref. 7, we adopt the symmetric-gate operation, in 
which the electric reference point of the new gate voltage Uq is set at the middle of the source- 
drain potential, Uq = Vq — Vd/2. Even if the external field is not applied, Vd = Uq = 0, 
electrons are so transferred from the low work-function side to the high work-function side 
that the chemical potentials coincide to reach equilibrium. There appears a scalar potential 
Vi from the redistributed charge density. Note that the scalar potential in electromagnetism 
is given by dividing it by — e. In order to use the periodic boundary condition for finite Vd, 
we introduce a time-dependent vector potential A{t) along the channel. Then the transfer 
integral acquires the Peierls phase, 

t iii+ i(t) = t* +hi (t) = j i i>i+1 (t) | exp [ij-A(t)\ , (2) 

where t denotes time, e the absolute value of the electronic charge, a the lattice constant, and 
c the light velocity. In order for the Poisson equation to contain only the scalar potential, V{, 
we use the space-independent vector potential, 

A(t) = -cJ dt'E D , (3) 

for the space-independent component of the electric field, 

E D = -V B /L , (4) 

with L being the number of sites. The gauge invariance of the motion of electrons is proved 
in Appendix. 

The scalar potential V{ obeys the Poisson equation on the discrete lattice, 

v i+1 - 2vi + Vi-i = -Vn ({rii) - n B i) , (5) 

where the parameter Vpj comes from the long-range Coulomb interaction and n-Qi from the 
background charge. Note that the Poisson equation itself is for three dimensions, but the scalar 
potential is assumed to depend on the position along the channel. If we take an orthorhombic 
system with lattice constants a, b, and c (not the light velocity here) and dielectric permittivity 
e (not the site energy here), the Gauss law gives Vpi = 47ra 2 /(6c) x e 2 /(47rea), which is 
proportional to the so-called nearest-neighbor repulsion strength and inversely proportional 



4/21 



J. Phys. Soc. Jpn. Full Paper 



(a) 








Gate insulator 




VJ2 


, -VJ2 










(b) 
W2 



mm 



-VJ2\ 



Fig. 1. (a) Device geometry, with source and drain electrodes on the left and right, respectively, and 
a gate electrode on top. (b) One-dimensional model, where the effect of the gate bias is included 
in a boundary value of the Poisson equation. 



to the dielectric permittivity e. The parameter Vpi is set at Vp c in the crystal and at Vp e in 
the electrodes, with Vp e <C Vp c . In the ideally metallic electrodes, Vp e would be zero. We use 
a small but finite value for it to stabilize numerical calculations. For the background charge, 
the neutrality condition, 

53<"i> =^2n Bi , (6) 

i i 

is crucial for the boundary conditions at the crystal-electrode interfaces to be compatible with 
those far from the interfaces. For Vq=0, the parameter n^i is set at n Bc uniformly inside the 
crystal and at n-Qe uniformly in the electrodes. They depend on the gate bias Uq and the drain 
voltage Vb- Concerning the site dependence of nsi for Vd ^ 0, we have compared the results 
with different assumptions: i) n-Qi in the crystal is a constant in the left half and another 
constant in the right half; and ii) n^i in the crystal is a linear function of i. In both cases, 
the site dependence of n-Qi in the electrodes is also set like that in the crystal. With either of 
these assumptions, n-Qi is numerically obtained with a solution Vi to the Poisson equation by 
imposing the self-consistency. These assumptions are found to produce quantitatively similar 
results. To show the data below, we use the assumption i). 

In principle, the gate bias Uq directly determines the potential of the gate electrode, 
not of the channel. The potential distribution in the gate insulator is determined by the 
Laplace equation, with boundary conditions set by the electrodes, which are coupled with the 
Poisson equation for the channel [Fig. 1(a)]. For a solution to the Laplace equation for the gate 
insulator, see for example ref. 2. The gate bias Uq consequently modulates the potential along 
the channel. For small | Uq |, the modulation of the potential at the midpoint would be linear 
with respect to Uq. Even if the modulation contains higher-order terms, they would not alter 
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the conclusion below on whether the field-effect characteristics are ambipolar or unipolar. 
Then we phenomenologically introduce the gate efficiency factor a 11 [Fig. 1(b)] instead of 
considering the transverse field. The details of the gate insulator including its shape and size 
are thus ignored. When the voltage at the gate electrode is varied in some range, the amount 
by which the potential changes at the midpoint of the channel is smaller than this range. 
Thus a is generally smaller than unity. For the thinner dielectric film (i.e., for the larger gate 
capacitance), the smaller gate-bias modulation causes a given amount of the drain-current 
modulation. 3 ' 4 This can be interpreted as due to the larger a. 

In such circumstances as described above, the gate bias Uq appears with the gate efficiency 
factor a in a boundary value of the Poisson equation, 

vl/2 -vq = -aU G + (f>, (7) 

where % = L/2 is set at the middle of the crystal and % = in the electrode is the furthest point 
from the crystal-electrode interfaces. Because Uq appears only here as a boundary value, we 
rewrite oUq as Uq for simplicity. Then the boundary condition is given by 

vl/2 ~vo = -Uq + <f> . (8) 

When the vector potential is added to the scalar potential, the total potential is given by 
Vi — Vd(i/L — 1/2), which becomes +Vd/2 deep inside the source electrode (i = 0), — Uq + <j) 
at the middle of the crystal (i = L/2), and — lb/2 deep inside the drain electrode (i = L) if 
v o = V L is arbitrarily set at zero. 

For each value of Uq, we first set Vb at zero and iteratively solve the eigenvalue equation 
for the ground state in the unrestricted Hartree-Fock approximation simultaneously with 
the Poisson equation for the potential distribution. After setting Vb at a finite value, we 
solve the time-dependent Schrodinger equation in the same approximation by decomposing 
the exponential operator, as performed in ref. 12. The self-consistency is imposed at each 
site and time both in the time-dependent Schrodinger equation and in the Poisson equation. 
Thus the potential distribution is obtained simultaneously with the current density in the 
non-equilibrium condition. 

The drain current Id is defined in the same manner as in ref. 13, 

/d = ~k C t(it{m ' (9) 

where the current density at time t, given by 

J W = ~\ [ it i+i.i( t ) c i+i,<r c i,<r ~ iti,i+i(t)cl a Ci+ lia , (10) 

i,cr 

is averaged over the period, < t < At with At = 27rL/(4Vb)- Its sign is set so that generally 
Id > (<)0 for Vb > (<)0. The time-evolution of the current density was compared with the 
exact one 13 for the pure Hubbard model (without electrodes) at and near half filling in a 



6/21 



J. Phys. Soc. Jpn. Full Paper 

wide range of Vb with the periodic boundary condition. The results are quite similar, sharing 
the qualitative characters. It is like teeth of a saw as a function of time for small Vb, while 
it is a sinusoidal function with a large amplitude for large Vb- Later we compare the drain 
current obtained by the Lanczos method and that obtained in the unrestricted Hartree-Fock 
approximation for small systems attached to electrodes. Both results show ambipolar field- 
effect characteristics in Mott insulators, but unipolar ones in band insulators. 

3. Results 

In the quasi-one-dimensional Mott insulator (BEDT-TTF)(F 2 TCNQ) used for the field- 
effect transistors, the intermolecular overlap along the channel is estimated at 0.019, 8 which 
corresponds to a transfer integral of about 0.2eV, while the polarized reflectivity spectra show a 
charge gap of about 0.7eV. 14 The work-function difference is estimated at 4> ~ — leV with gold 
electrodes. 7 For numerical calculations, we use t c =l, as a unit, and U=2 for Mott insulators, 
which produces a charge gap of about Acg=0.68 at half filling in the present approximation. 
In most of the calculations below, we use t e =l. Thus, compared with the experimental values, 
the transfer integral in the crystal is larger, while that in the electrodes is smaller. We use them 
simply because the numerical convergence is generally stable and rapid for comparable transfer 
integrals. As shown later, the main conclusion is not altered by different transfer integrals. 
Below we frequently use Vp c =0.05 and Vp e =10 -3 , so that Vp c /Vp e =50. The numerical results 
are insensitive to Vp c or Vp c if the ratio is unchanged, so that we mention the ratio only. It 
is shown later that different ratios do not alter the main conclusion. Below we mostly use the 
100-site periodic system with the central 51 sites belonging to the crystal and the peripheral 
49 sites to the electrodes. 

When the work function of the crystal and that of the electrodes coincide (<f)=0), the drain- 
current (Id )-gate- voltage (Ug) characteristics are symmetric with respect to the exchange of 
Uq and —Ug owing to the electron-hole symmetry. This property is irrespective of whether 
the electronic state is a Mott insulator [Fig. 2(a)] or a band insulator [Fig. 3(a)]. 

Mott insulators show nearly symmetric Id-Ug characteristics for finite work-function dif- 
ferences (Fig. 2). For small Vb and | Uq |< Acg/2, In is suppressed by the charge gap. 
The drain current Id increases with electron injections (Ug > 0) and with hole injections 
(Ug < 0) in a similar manner for any Vb- For large Vd, further carrier injections (with in- 
creasing | Ug I) lead In to a maximum and slight reduction after reaching the maximum. This 
is due to the small bandwidth of the electrodes used here. In fact, this behavior disappears for 
electrodes with a wide bandwidth, as shown later. It is not reported in the experiment with 
gold electrodes as well. 7 

Around Ug=0, a threshold exists in Vb above which Id becomes finite and rapidly increases 
with Vb for Mott insulators. The threshold Vbth is in the range 0.8 < Vbth < 0.9 for 4> = 
[Fig. 2(a)], 0.7 < Vb th < 0.8 for cp = -0.5 [Fig. 2(b)], 0.6 < F Dth < 0.7 for = -1.0 



7/21 



J. Phys. Soc. Jpn. 




Full Paper 




Fig. 2. I-d-Uq characteristics at various drain voltage Vb for Mott insulators with U=2 and 8t=0. 
The work-function differences are, (a) </>=0, (b) <j>~— 0.5, (c) <f)——1.0, and (d) (/>=— 1.5. The other 
parameters are t c =t c =l, Vp c /Vp e =50, and L=100. 

and -1.5 [Figs. 2(c) and 2(d)], and 0.5 < Vbth < 0.6 for (f> = -2.0 (not shown). The 120- 
site system with 63 crystal sites and 57 electrode sites has slightly different thresholds, e.g., 
0.7 < Vbth < 0.8 for (f) = —1.0 (not shown), but the tendency is unchanged. The threshold 
Vbth moderately decreases with increasing work- function difference. Above the threshold Vbth) 
the I-q-Uq characteristics generally deviate from the symmetric ones largely compared with 
those below Vbth- There appears no general rule concerning the asymmetry (e.g., which of the 
electron injections and the hole injections are more effective in increasing Ir>) above Vbth- The 
asymmetry depends on the bandwidth of the electrodes, as shown later. The largest In value 
at the optimum Uq at a given Vb > Vbth moderately decreases with increasing work-function 
difference (see the vertical scales). In general, the larger work-function differences bring about 
the energetically higher and spatially wider Schottky barriers, so that the charge gap is more 
easily collapsed by the drain voltage (thus Vbth is reduced) and the drain current is generally 
reduced. 
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Below Vbth; however, we find a general tendency (as shown here and below with different 
parameters) concerning the deviations from the symmetric I-q-Uq characteristics. When the 
work function of the crystal is lower than that of the electrodes {(ft < 0) , Id is slightly larger 
for Uq < than for Uq > for a given | Uq | at Vb < Vbth, as if the hole mobility is 
higher than the electron mobility. This is consistent with the experimental finding. 7 It can 
be concluded from this result with the electron-hole exchange operation that, for (ft > 0, 
the I-q-Uq characteristics behave as if the electron mobility is higher. The drain current Id 
sometimes shows an oscillating dependence on Uq especially for small Vb- It is caused partly 
by averaging the current density, which generally oscillates with a C/c-dependent frequency, 13 
over the fixed time period, and partly by iteratively obtaining the potential distribution during 
the application of the drain voltage. It does not have any physical meaning. 

The drain current Id has a minimum at Uq = Uq 111 . The minimum point U™ m is exactly 
zero for (j)=0, but it generally deviates from zero. The dependence of Uq 1u on (ft is not mono- 
tonic. It is shifted to a positive value for small | (ft | and then to a negative value for large | (ft \ 
in Fig. 2 for (ft < 0. The sign of Uq" 1 can be altered by using a different bandwidth of the 
electrodes, as shown later. Thus the behavior of Uq 111 is not simple. In the experiment, Uq 1u 
is reported to be 2-3V. 7 If the gate efficiency factor were unity, it would be much smaller. 

We have a few comments on the experimental observation. The nonlinearity in Ib(Uq ^ 
0) — Id(Uq = 0), as a function of Vb, becomes more prominent at low temperatures. 7 The 
current-voltage characteristics deviate from those expected in the gradual channel approxima- 
tion largely there. In fact, the assumption on which this approximation is based does not hold 
in the material with the quasi-one-dimensional and coherent band transport. This point will 
be discussed in the last section. At high temperatures, however, thermal fluctuations obscure 
the intrinsic nonlinearity and would make the characteristics similar to those expected in this 
approximation . 

In order to compare the field-effect characteristics in Mott insulators with those in band 
insulators, we used the tight-binding model with alternating transfer integrals, that with 
alternating site energies, and that with staggered magnetic fields. We show numerical results 
for the model with alternating transfer integrals only because the results for the other models 
are very similar. We use <5i=0.17, which also gives a charge gap of Acg=0.68 at half filling. 

Band insulators clearly show very asymmetric I-q-Uq characteristics for finite work- 
function differences [Figs. 3(b)- 3(d)]. Only when the work functions coincide, Id is suppressed 
only for | Uq |< Acg/2 and suddenly increases with Uq or Vb for | Uq |> Acg/2 [Fig. 3(a)]. 
As shown later in Fig. 4, the Schottky barriers with (ft < are higher for the electron injections 
Uq > than for the hole injections Uq < at a given | C/q I - The drain current Id increases 
with the hole injections, but it remains suppressed for the electron injections. Although Id 
increases with Vb, this behavior is unchanged. A threshold in Vb might be defined as long as 
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Fig. 3. Id-Uq characteristics at various drain voltage Vb for band insulators with U=0 and <5i=0. 17. 
The work-function differences are, (a) </>=0, (b) 4>=— 0.5, (c) (j>=— 1.0, and (d) 4>—— 1.5. The other 
parameters are the same as in Fig. 2. 



we focus on Id near [/q 111 only, but it is insignificant. These unipolar field-effect characteristics 
in band insulators are generally obtained irrespective of whether the charge gap is introduced 
by alternating transfer integrals, alternating site energies, or staggered magnetic fields. They 
are in marked contrast to the ambipolar field-effect characteristics in Mott insulators, where 
the charge gap is introduced by the electron correlation. 

The asymmetric characteristics in band insulators with various work-function differences 
(4> ^ 0) are similar if the polarity of the gate bias is fixed ((f) < in the present case). 
With increasing the work- function difference | (j> \i Id f° r the gate-bias polarity with higher 
Schottky barriers (Uq > in the present case) is further suppressed, while Id with lower 
Schottky barriers (Uq < in the present case) is further enhanced [Figs. 3(b)-3(d), see the 
vertical scales]. 

The spatial distribution of the total potential vi — Vo(i/L — 1/2) is shown in Fig. 4(a) 
for the Mott insulator, and in Fig. 4(b) for the band insulator. Because the work function 
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Fig. 4. Potential, Vi — Vh(i/L — 1/2), as a function of the site index, i, at positive and negative gate 
biases Uq, (a) for the Mott insulator with U=2 and St=0, and (b) for the band insulator with 
£7=0 and 5t=0.17. The crystal sites are 25 < i < 75. The other parameters are 4>= — 1.0, Vb=0.6, 
t e =t c =l, and Vp c /Vp o =50. 

of the crystal is set lower than that of the electrodes, electrons transfer from the crystal 
to the electrodes, which make the second derivative of the potential positive in the crystal 
and negative in the electrodes, lowering the Fermi level of the crystal relative to that of the 
electrodes. For positive gate bias, Uq > 0, the Fermi level of the crystal is further shifted 
downwards to accommodate electrons, so that the Schottky barriers at the interfaces are 
made higher. Meanwhile, for Uq < 0, the potential is so modified as to accommodate holes, 
which make the Schottky barriers lower. In both cases, the potential has a long tail toward 
the middle of the crystal, which is in contrast to the MOSFET, where the larger impurity 
concentration shortens the depletion length. Roughly speaking, the potential distribution in 
the Mott insulator is similar to that in the band insulator. However, if one looks carefully at the 
potential difference between positive and negative gate biases, that for the Mott insulator is 
found different from that for the band insulator, especially near the crystal-electrode interfaces. 
In the case without gate bias (Ug=0), the spatial variation of the density of states in the 
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Fig. 5. Density of states in the crystal, 25 < i < 75, without gate bias, Uq=0, (a) for the Mott 
insulator with U=2 and 6t=0, and (b) for the band insulator with U=0 and 5t=0.17. The other 
parameters are the same as in Fig. 4. 

crystal is shown in Fig. 5(a) for the Mott insulator, and in Fig. 5(b) for the band insulator. In 
both cases, the gap is formed, so that the drain current does not flow. The upper and lower 
bands correspond to the upper and lower Hubbard bands in Fig. 5(a), and to the conduction 
and valence bands in Fig. 5(b). At the bottom of the upper bands and at the top of the lower 
bands, the density of states increases abruptly, forming ridges. The ridges have a curvature 
according to the spatial variation of the potential shown in Fig. 4. The space and energy 
dependence of the density of states in the Mott insulator is also similar to that in the band 
insulator as long as the drain current does not flow. Because their static properties are similar, 
the difference in their field-effect characteristics must be understood from the dynamical point 
of view. 

In the case with a finite gate bias (C/g=0.2), the spatial variation of the density of states in 
the crystal is shown in Fig. 6(a) for the Mott insulator, and in Fig. 6(b) for the band insulator. 
Here, the gap is formed in the band insulator only. The parameters used in Figs. 6(a) and 6(b) 
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Fig. 6. Density of states in the crystal, 25 < i < 75, with a gate bias, Z7g=0-2, (a) for the Mott 
insulator with U=2 and 6t=0, and (b) for the band insulator with U=0 and 5t=0.17. The other 
parameters are the same as in Fig. 4. 

here are the same as those in Figs. 2(c) and 3(c), respectively, with Vb=0.6. The drain current 
flows in the Mott insulator here. Precisely speaking, the Hubbard system with the finite gate 
bias and the large-enough drain voltage is conducting. Thus the density of states for the Mott 
insulator becomes rather flat and featureless. The dynamic properties are determined in a 
self-consistent manner from the electron correlation and the Schottky barriers. This point will 
be discussed later. 

In order to see how universal the field-effect characteristics of the Mott insulators are, we 
used different parameters. They have turned out to be always ambipolar. In Fig. 7(a), the 
bandwidth of the electrodes is set larger than that of the crystal. Compared with the figures 
already shown, the drain current shows rather irregular behavior especially for large Vb- This 
is due to numerically slower and poorer convergence generally for large bandwidths of the 
electrodes. Nevertheless, it is clear that the maximum of Id as a function of Uq for large Vb 
in Fig. 2(c) disappears in the present case from the Uq range shown here. In other words, the 



13/21 



J. Phys. Soc. Jpn. 



Full Paper 




Fig. 7. I-d-Uq characteristics at various drain voltage Vb for Mott insulators with U=2 and 8t=0. 
The bandwidth ratios are, (a) t e /t c =2, and (b) t c /t c — 0.5. The other parameters are </>=— 1.0, t c =l, 
V Pc /V Po =50, and L=100. 

density of states in the electrodes substantially influences the behavior of Id as a function of Uq 
for large Vb- Compared with Fig. 2(c) for t c /t c =l, the drain current is generally reduced by the 
larger t e here. It indicates that the mismatch of the bandwidths would contribute to backward 
scatterings of electrons entering the crystal, reducing the current flowing in the crystal. In 
Fig. 7(b), the bandwidth of the electrodes is set smaller than that of the crystal. Here again 
the drain current is reduced, compared with the case of t e /t c =l. Obvious differences between 
Figs. 7(a) and 7 (b) are found in the deviations from the symmetric Id-Uq characteristics 
above Vbth- I n the former the hole injections are more effective in increasing Id, while in the 
latter the electron injections are more effective. The minimum points C7™ in around Vbth are 
also shifted in the opposite directions from the point at t c /t c =l. 

In Fig. 8, the ratio of the Coulomb parameters, Vp c /Vp e , is varied. Because of numerically 
poorer convergence, we use a smaller system than before, but the drain current shows rather 
irregular behavior. With increasing the ratio Vp c /Vp e , the I-q-Uq characteristics appear to 
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Fig. 8. Id~Uq characteristics at various drain voltage Vb for Mott insulators with U=2 and 6t=0. 
Different ratios of the Coulomb parameter in the crystal to that in the electrodes, Vp c /Vp e , are 
used for comparison. The other parameters are <f)=—l, t e =t c =l with smaller L=80 for convergence 
reasons. 



approach symmetric ones. In the realistic situation, this ratio would be much larger. In any 
case, the ambipolar field-effect characteristics are robust for the Mott insulators. We have 
shown the results with U=2 for the Mott insulators. We checked the robustness by using 
different strengths of on-site repulsion. For instance, a weak on-site repulsion leads to a smaller 
charge gap, so that all the relevant energies including Uq and Vb are scaled down. The It>-Uq 
characteristics are always nearly symmetric below Vbthj demonstrating the robustness. 

In order to clarify the origin of the robustness of the ambipolar field-effect characteristics, 
we show snapshots of charge density (n») for both gate bias polarities in Fig. 9 before the drain 
voltage is applied. It is clearly seen that electrons transfer from the interfacial region inside 
the crystal to the whole region of the electrodes so as to adjust the relative Fermi levels when 
the work function of the crystal is set lower than that of the electrodes. Near the middle of the 
crystal, the charge density is almost flat and near unity due to the strong on-site repulsion. 

For Uq < the Schottky barriers are lower [Fig. 4(a), upper curve], and the depletion 
length is longer [Fig. 9(a)], while for Uq > the Schottky barriers are higher [Fig. 4(a), 
lower curve] and the depletion length is shorter [Fig. 9(b)]. For the gate-bias polarity with 
the higher Schottky barriers {Uq > 0), the deviation of the charge density well inside the 
insulator from unity is always larger than that for the opposite gate-bias polarity, indicating 
that the umklapp scattering is weaker. Therefore, the robustness of the ambipolar field-effect 
characteristics is a consequence of balancing the correlation effect with the barrier effect in a 
wide parameter space. 

So far, we have used the unrestricted Hartree-Fock approximation to replace the many- 
body wave function by the self-consistent Slater determinant of the one-body wave functions. 
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Fig. 9. Snapshots of charge density, (a) for Uq~— 0.2, and (b) for [7 G =0.2. The sites 25 < i < 75 
are for the Mott insulator with U=2 and 5t~0. The other parameters are (/>=— 1.0, t =i c =l, and 
W^p o =50. 



In order to confirm that the ambipolar field-effect characteristics of the Mott insulators are 
not artifacts of this approximation, we employ the Lanczos method to treat the many-body 
wave function exactly. The potential distribution is obtained first in this approximation and 
averaged over the period, < t < At with At = 2ttL/(4Vd). With this static potential 
distribution, the time-dependent Schrodinger equation is solved by the Lanczos method. In 
the 10-site periodic system with 7 sites belonging to the crystal and the rest to the electrodes, 
5 up-spin and 5 down-spin electrons are included. In this small system, a band insulator 
with a small gap shows rather symmetric I-d-Uq characteristics if the localization length 
is comparable with or larger than the system size. Then, we compare the Mott insulator 
with U=4 and the band insulator with 8t=0.77, both of which have a large charge gap of 
about Acg=3.1, and use the work-function difference <p=—2. The band insulator shows very 
asymmetric Ib-Uq characteristics again: Id for Uq=— 2 is about four times as much as Id 
for Uq=2 when Vb=2 is applied (not shown). On the other hand, the Mott insulator shows 
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Fig. 10. Id-Ug characteristics for Mott insulators with U=A and St=0, (a) calculated with the help 
of the Lanczos method, and (b) calculated in the unrestricted Hartree-Fock approximation. The 
parameters are 4>~— 2.0, Vb=2.0, t c =t c =l, Vp c /Vp e =50, and L=10. 



nearly symmetric Ib-Uq characteristics, both by the Lanczos method [Fig. 10(a)] and in the 
Hartree-Fock approximation [Fig. 10(b)]. Therefore, the ambipolar field-effect characteristics 
are the intrinsic properties of Mott insulators. The magnitude of the charge gap is about 
Acg=3-1, both by the exact diagonalization and in the present approximation. The origin of 
the difference between the Uq dependence of the exact Ijy and that of the approximate one is 
unclear. 

4. Summary and Discussion 

In order to clarify the mechanism of the ambipolar field-effect characteristics in MISFET 
device structures based on organic single crystals of the quasi-one-dimensional Mott insulator 
(BEDT-TTF)(F2TCNQ), 7 we have performed model calculations taking account of the fact 
that the gate bias does not only modulate the carrier density but also affect the Schottky bar- 
riers in this low-dimensional system showing coherent band transport. 7 The one-dimensional 
Hubbard model is used for Mott insulators and attached to the tight-binding model for source 
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and drain electrodes with a generally different work function. The scalar and vector potentials 
satisfying the Poisson equation are added to the model, which originate from the long-range 
Coulomb interaction and form the Schottky barriers. The drain voltage and the gate bias mod- 
ulate the boundary values of the Poisson equation. The time-dependent Schrodinger equation 
is solved in the unrestricted Hartree-Fock approximation, simultaneously with the Poisson 
equation. The self-consistency is imposed at each site and time in both equations. 

Mott insulators show ambipolar field-effect characteristics in a very wide parameter space 
spanned by the work-function difference, the bandwidth difference, and the Coulomb param- 
eters in the Poisson equation. Thus, the experimentally observed, ambipolar carrier injections 
are indeed general characteristics of quasi-one-dimensional Mott insulators. For the gate-bias 
polarity with higher Schottky barriers, the correlation effect is weakened accordingly. In other 
words, the robust ambipolar characteristics are caused by balancing the correlation effect with 
the barrier effect. This is in marked contrast to the fact that band insulators generally show 
unipolar field-effect characteristics as long as the work function of the crystal is different from 
that of the electrodes. 

The difference in the transport property can be understood as due to the difference in the 
collectivity. To contribute to the current density, carriers must go through the crystal-electrode 
interfaces. When the crystal is a band insulator, each carrier is affected by the Schottky barrier 
independently of other carriers. As the barrier becomes higher, it becomes more difficult for 
the carrier to go through it, so that the current density is reduced. When the crystal is a Mott 
insulator, all the electrons are correlated, repelling nearby electrons. If the electron density 
deviates from the equilibrium one in such a way that the positive or negative charge density is 
locally accumulated, it is energetically unfavorable. Thus the electrons move collectively. This 
is why the transport property is insensitive to the local variation of the potential distribution 
near the interfaces in Mott insulators. The charge density distribution is determined in a 
self-consistent manner to satisfy both the Schrodinger equation and the Poisson equation. 
Thus the influence of the on-site repulsion and that of the long-range Coulomb interaction 
counterbalance each other. Balancing the correlation effect with the barrier effect causes the 
robust ambipolar characteristics in Mott insulators. 

Such collective transport is reminiscent of the knock-on mechanism in ion channels. 9 ' 10 
K + ions hop in single file (i.e., in a queue) from one binding site to the next as permeation 
proceeds. Although hopping processes occur temporally in a random fashion, they are spatially 
correlated. The neighboring ions hop together. Roughly speaking, the occupation of ions in the 
binding sites is either 0101 or 1010, which is reminiscent of a charge-ordered Mott insulator. 
The effective electrostatic repulsion between the ions is manifested at short distances. The 
concerted multi-ion conduction mechanism is determined largely by the multi-ion free energy 
surface rather than the dissipative and frictional forces. Thus the ion translocation in the 
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multi-ion potential energy surface is essential for achieving a high throughput rate, as clarified 
by a modified Poisson-Boltzmann theory. 9 ' 10 

When the motion of carriers is coherent and restricted to a one-dimensional system, the 
gradual channel approximation does not generally work because the contact resistance domi- 
nates the channel conductance. Thus the barrier effect is relatively important. In this sense, 
the mobility reported in the experiment on the quasi-one-dimensional Mott insulator, 7 is an 
effective one because the mobility is estimated by fitting the current-voltage relation to that 
based on the gradual channel approximation. The fact that the mobility increases with the 
drain voltage and depends on the polarity of the gate bias indicates that it is an effective 
one. The estimated mobility values are indeed very low (by some three orders of magnitude), 
although the coherent band transport is expected from the temperature dependence of the 
mobility. 7 The coherent band transport would enhance the relative importance of the Schottky 
barriers. 

It is true that the intermolecular overlaps along the molecular long axis are finite, although 
they are very small. Considering the transverse field (i.e., the transverse distribution of the 
potential) and the consequent transverse motion of electrons are necessary not only when a 
crossover from one-dimensional to three-dimensional systems is discussed but also when the 
gate bias is so large that the modulation of the longitudinal potential distribution, especially 
near the crystal- (source/drain) electrode interfaces, is not described by the gate efficiency 
factor. 
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Appendix: Gauge Invariance of Electron Dynamics 

For simplicity, we treat the one-dimensional model containing the scalar-potential and 
transfer terms only, by setting e% = Ui = and | tj |= to i n ec L- (!)• The total force from 

the scalar and vector potentials is shown to be given symbolically by 

-e.Etot = -eE D - d x v , (A-l) 

independently of the gauge. For intuitive understanding, let us first take a plane wave, exp[i/cr] 
at the position r = ja with integer j. The Peierls phase, exp [ieaA(t) / (he)} = exp [—ieEnat/H], 
shifts its momentum hk by —eE^t, which corresponds to the force —eE-£> from the vector 
potential. Thus the above equation looks reasonable. 
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To prove the above equation, we consider how the expectation value of a (generally off- 
diagonal) one-body density, (exp (ineEj^at/h) cj +ra a Cj iCr ), evolves in time. It is straightforward 
to derive the equation, 

ihd t (ex.p [ineEvat/h] c] +njff c,- j(T ) 

= [-neE-oa - {v j+n - vj)} exp [meE^at/h] (Cj +n a c jjCr ) 

-to exp [i(n - l)eE D at/h] {c] +n (J c j+1 , c - c) +n _ l a c jja ) 

+t exp [i(n + l)eE D at/h] {c) +n+l a Cj^ - c) +n a Cj- h(T ) . (A-2) 

This shows that the time-derivative of (exp (ineE^at/h) ct +n a Cj j(7 ) for any integer n is always 
given by a linear combination of themselves for n = 0, 1, • ■ ■ , L— 1. All the coefficients are either 
±io or [ — neE^a — {vj +n — Vj)\. Thus their time-evolution is governed by the contribution from 
the vector potential, —neE^a, and that from the scalar potential, — (vj +n — Vj), which appear 
always as the sum of them. Symbolically, it corresponds to na times the quantity —eE-^ — d x v 
in the continuum limit. Thus the above equation is proved. 

This proof is easily extended to the case, in which the vector potential depends 
on the bond on which a finite transfer integral is assigned. The time-derivative of 
(exp (ie(E Dj+n _ 1/2 H h E Bj+1/2 )at/ti) c^ j+n ^c ji(T ) for any integer n is again given by a lin- 
ear combination of themselves for n = 0, 1, • • • , L — 1. All the coefficients are then either ±io 
or [— e(E- D j +n _i/2 + • • • + Ej)j + i/2)a — (vj +n — Vj)~\ . If we take a space-dependent vector po- 
tential, it modifies the Poisson equation. Thus the convenient gauge is the space-independent 
vector potential that absorbs the finite Vb so that the scalar potential is periodic and the 
vector potential is not involved in the Poisson equation. 
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